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1 Introduction 

Fluid dynamics [1, 2] is an effective long-wavelength description of most classical or quan¬ 
tum many-body systems at nonzero temperature. For neutral fluids in flat space, the 
hydrodynamic equations are derivable from conservation of the fluid’s stress-energy tensor 



( 1 . 1 ) 


( 1 . 2 ) 


where e, u M are the fluid’s energy density and four-velocity held, whereas rj llh , stands for 
Minkowski metric tensor. The pressure P is specified through equation of state P = P(e), 
calculable from underlying microscopic theory. Deviations from thermal equilibrium are 
collectively encoded in dissipation tensor 



(1.3) 


where = r + u^u u is a projector on spatial directions. 

At each order in derivative expansion, 11^ is fixed by thermodynamics and symme¬ 
tries, up to some transport coefficients. The latter have to be calculated from microscopic 
description of the fluid rather than from hydrodynamics itself. In what follows, we focus 
on conformal fluids in 4D Minkowski spacetime, so the condition = 0 implies £ = 3 P. 
The first order derivative expansion gives the Navier-Stokes term 


iv = -2Vo d P u - 


'0 u H Uj V’> 


(1.4) 
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where r/ 0 is a shear viscosity. At second order, there are five additional transport coeffi¬ 
cients [3, 4]. 

AdS/CFT correspondence [5] relates strong coupling physics of gauge theories with 
large number of colors N to weakly coupled gravity in (asymptotic) AdS space. As a partic¬ 
ular example, it maps hydrodynamic fluctuations of a boundary fluid into long-wavelength 
gravitational perturbations of a stationary black brane in asymptotic AdS space [6-8]. Vis¬ 
cosity and all other transport coefficients could be computed from the gravity side of the 
correspondence. The ratio of rj 0 over the entropy density s was computed in [6, 7, 9] 


r h = J_ 

S 4-7T 


(1.5) 


and was found to be universal for all gauge theories with Einstein gravity duals [10-12], 
The value (1.5) was further conjectured to be Nature’s lower bound for y 0 /s [13]. 

The relativistic Navier-Stokes hydrodynamics is well known to violate causality, that is 
it admits propagation of signal faster than the speed of light. Inclusion of any finite number 
of additional derivative terms in would not render the theory into causal. All-order 
derivative resummation is necessary to restore causality. In [14-16], we built upon the work 
of [17] and linearly resummed derivative terms (see [18-22] for boost invariant case) for 
fluids dual to pure Einstein gravity. In a parametrically controllable approximation, where 
we only collect terms linear in amplitude of the fluid velocity, has a compact form, 


IV = —277 [u a d a , d a d a \ 


c [ u a d a , d a d a } d^d v d a u a . 


( 1 . 6 ) 


Here r] and £ are derivative operators, which upon expansion in a series would generate 
the usual gradient expansion. Thanks to linearization, we can study these operators in 
Fourier space, via replacement d^ —> (—ioj,iq ). Then the operators rj and C, are turned 
into functions of momenta and are referred to as viscosity functions. In momentum space, 
the constitutive relation (1.6) is 

Tl^v(u, q) = -2 r)(u, q 2 )iq^u u (u, q) + C(w, q 2 )iq fi q u q a u a (uj, q). (1.7) 


The viscosity functions r/ and Q were computed exactly in [14-16] and were observed to 
vanish at very large momenta, signaling restoration of causality in the dual CFT. For self- 
consistency of presentation we will flash these results in section 3 below. 

Vanishing of the viscosities at large frequencies is a necessary condition for causality 
restoration. To better understand the physical role of the viscosity functions, we turn them 
into memory functions via inverse Fourier transform of (1.7) 


n 




/ oo r i 

dt' 2 fj(t - t', q 2 )d fJ ,u u (t') + ((t - t', q 2 )d^d u d a u a (t') , (1.8) 

-OO L J 


where 


V(t,q 2 )= [ -^=ij(uj,q 2 )e lujt , C(t,g 2 )= / 4^C(^,g 2 )e 

J — oo V J — oo V 
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(1.9) 






Here fj(t,q 2 ) and ((t,q 2 ) are memory functions in mixed-(t, q 2 ) representation. Causality 1 
requires memory functions to have no support in the future: fj(t — t') ~ 0(t — t') and 
C (t-t) r\ j 0(t — t'). In other words, the current n ^(t) at time t should be affected by the 
state of the system in the past only. So, for a causal theory becomes 


n 



t',q 2 )d^Uv(t') + C(t 


t',q 2 )d^d u d a u a (t’) 


( 1 . 10 ) 


As has been discussed in [16], a typical memory function-based formalism [2, 24], as a 
phenomenological model, would set the low limit of integration in (1.10) to zero, turning 
thus defined hydrodynamics into a well-posed initial value problem. 

In [16], the memory function fj[t) was evaluated from exact computations in the dual 
Einstein gravity. It was indeed found to be proportional (up to numerical noise) to 0(i) as 
could be seen from Figure 1 (in units ttT = 1). 



Figure 1. Memory function fj(t,q 2 ) evaluated in [16] for hydrodynamics dual to pure Einstein 
gravity. Left: 3D plot as function of t and q 2 . Right: 2D plots as functions of time t: different 
curves display results with different q 2 (from the rightmost: q 2 = 0,1, 2, 3). 

Beyond N —> oo and ’t Hooft coupling A —> oo limits, the ratio (1.5) gets corrected. 
Finite N or A corrections arising from stringy or quantum effects introduce, beyond Einstein 
gravity, terms with higher derivatives of curvature. Exact forms of these terms generated in 
string theory are not known in general: the first higher derivative correction is expected to 
be the curvature squared. Of particular interest is a ghost-free Gauss-Bonnet combination, 
which generates equations of motion of second order only. Adding a Gauss-Bonnet term to 
the gravitational action is equivalent to introducing some 0(1/N) corrections in the dual 
gauge theory, whereas the Gauss-Bonnet coupling a is related to the difference between two 
central charges of the dual CFT. 

From the string theory point of view, the Einstein-Gauss-Bonnet (EGB) gravity should 
be considered as phenomenological effective low energy theory. One may, however, con¬ 
sider EGB gravity on its own, as a UV complete theory. Still applying the rules of the 
AdS/CFT correspondence, one finds that the Gauss-Bonnet correction violates the lower 

1 In the limit IV —> oo and A —> oo, causality of M — 4 super-Yang-Mills plasma was analyzed [23] by 

studying pole structures of retarded correlators. 
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bound (1.5) [25, 26] (see also [27-32]). Non-perturbative Gauss-Bonnet corrections to sec¬ 
ond order transport coefficients in conformal fluids dual to EGB gravity were considered 
in [33-35]. Furthermore, causality of the dual CFT sets constraints on possible values of the 
Gauss-Bonnet coupling. In [26, 36-38], the coupling a (stripped of units) was constrained 
to be 

7 u 

( 1 . 11 ) 


-< a < 

72 ~ ~ 


9 

200 ’ 


where the lower (upper) bound was obtained by requiring the front-velocity in the sound 
(scalar) channel of the dual CFT not to exceed the speed of light. Positivity of energy flux 
in thought experiments done in conformal colliders [39] also constrains values of a [40-42], 
Remarkably, constraints on a from causality and positivity of energy flux were found to 
match [38, 40-42], Stability of the dual plasma also sets constraints on a. [43, 44], More 
recently, causality violating effects due to higher derivative corrections to Einstein-Hilbert 
action were discovered in high energy scattering processes of gravitons off shock waves [45] 
and strings off branes [46]. Pure EGB gravity was concluded to be a-causal for a of order 
one. Causality is restored by adding an infinite tower of extra massive particles with spins 
higher than two [45]. 

In this work, we would like to explore the effects of the Gauss-Bonnet corrections on 
transport coefficients, beyond known results at first and second order. To this goal we 
consider hydrodynamics dual to EGB gravity and calculate Gauss-Bonnet correction to 
viscosity functions rj and £. Given previous constraints on a, we limit our study to small a 
only. To linear order in a, the fluid’s energy density and pressure are 


e = 3P = 3(1 + 3a) (vrT) 4 . (1.12) 

The entropy density is evaluated from s = dP/dT 

s = 47r(l + 3a) (vrT) 3 . (1-13) 


In the hydrodynamic limit, the viscosity functions are expandable in momenta, 


r](uj, q 2 ) = (1 — 5a) + - [(2 — In 2) — (21 — 5 In 2) a] ioo — 


— (67t — 7t 2 + 24 
48 v 


-36 In 2+ 12 In 2 2) - 6.53(280)a1 u? - ( - - 2.11(320)a) q 2 + • • • , (1-14) 


C(w, q 2 ) =— [(5 - vr - 2 In 2) + (15vr - 87 + 30In 2) a] + 


where the fluid’s temperature is normalized to ttT = 1 and all the momenta are set to be 
measured in these units. For a positive a, the first term in rj yields violation of the viscosity 
to entropy bound [25, 26, 47, 48], 

— = T-(l-8a). (1.15) 

S 47T 

The second term in r] is the relaxation time, calculated in [33, 49]. The remaining terms 
are new third order transport coefficients. The underlined terms are our numerical results 
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for ct-corrected pieces. To resum the derivative terms to all orders, we numerically compute 
the viscosity functions for generic u and q 2 . The viscosity functions are formally expanded 
in a 

Tj = r/°) + a 7 / 1 ) + O (ct 2 ) , C = C (0) + a C (1) + O (a 2 ) j (1.16) 

where rand are the viscosity functions computed for pure Einstein gravity in [14, 15]. 
The results of this calculation, particularly new results on rand (4 1 ), are presented 
in subsection 3.2.2. When Fourier transformed into memory functions, we find that the 
Gauss-Bonnet correction (and also (p 1 )^)) i g also vanishing at negative times (see 

Figure 9) 2 . 

In section 2, we present the holographic setup. A boosted black hole solution of the EGB 
gravity in asymptotic AdSs space is introduced. Following [4], gravitational perturbation is 
induced by locally varying boost velocity and black hole temperature. We then parameterize 
additional bulk metric corrections in terms of ten functions h, k, ji and ctjj, which are 
both functions of holographic coordinate and functionals of the fluid velocity u p . The 
boundary stress-energy tensor is read off from holographic renormalization, being expressed 
in terms of near-boundary behavior of h, k, ji and ctij. In section 3, we solve the Einstein 
equations for the metric corrections. Thanks to linearization in the velocity amplitude, all 
bulk metric corrections can be decomposed in the basis formed from rq. As a result, in 
Fourier space, the Einstein equations turn into second order ordinary differential equations 
for decomposition coefficients. Solutions to these equations reveal the information about 
the viscosities. We then discuss effects of the Gauss-Bonnet correction on the viscosity 
functions. Section 4 is devoted to summary and discussion. Some computational details 
are provided in Appendix A. 

2 Holographic setup for Einstein-Gauss-Bonnet gravity 

Our representation is largely based on [50]. We start from the EGB gravity with a negative 
cosmological constant A = —6/l 2 in 5D spacetime manifold At, 

S = n — f d 5 Xy/—g (i? — 2A + al 2 Lq b) + S sm: + S c , t., (2-1) 

1d7TG7V Jm 

where the Gauss-Bonnet term Lq b is 

Lqb = RmnpqR MNFQ ~ 4 RmnR MN + R 2 - (2.2) 

We use a mostly plus signature for the bulk metric gMN■ To have a well-defined variational 
principle, the surface term S sur computed in [51, 52] was added to (2.1), 

Ssur = TT-^ [ 7 [K - 2cd 2 (J + 2 G^K^)\ , (2.3) 

ovtGat Jqm 

where the first term is the Gibbons-Hawking surface action. The tensor J^ v is defined as 
•V = -3 (2 KK^KP + K pa K^K iW - 2K lip K^K a „ - K 2 K pv ) , (2.4) 

1 In the first version of this preprint, we made a wrong statement on causality violation based on numerical 

Fourier transform, which was later found to be lacking sufficient accuracy. 
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where K^ = 'y®'V a n u and 7 ^, n p are the induced metric/outing normal vector on/to a 
constant r-slice dM. The Einstein tensor G pv and are compatible with 7 ^. 

In asymptotic AdS space, UV divergences near conformal boundary can be removed 
by holographic renormalization [53, 54], For the EGB gravity, the counter-term action S c ,t. 
was first constructed in [50] following previous studies [55-57], 


Set. = 


1 


8 ttGn Jqm 

with the coefficients <5i and 62 having the forms 


d 4 xv' = 7 ($1 ~ yftpy]) , 


(2.5) 


— 1 — 8 a + y/\ — 8 a a ~>0 3 a , 2 . 

61 = 7- / —► - J + T + °( a )> 

V Aal 2 v 1 — y/1 — 8 a ^ ^ 

\/4a / 2 (3 — 8 a — 3\/l — 8 a) Q _> 0 l 31 /n . 2 . 

5 2 = --- n3/2 - > - + -a + 0(a ). 

2 (1 - Vl - 8 a) ^ 2 2 


( 2 . 6 ) 


Up to a conformal factor, the stress-energy tensor of the boundary CFT is obtained by 
varying (2.1) with respect to 7 /ll/ . The boundary stress-energy tensor is [50], 


T/iv = lim T tiv (r) 


= — lim ——— 
r —>00 8ttGn 


Kfiv K ltlu 5i7/ii/ ^ 2 Gfin 2al [ qQ'Yhv 


The tensor is dehned as 


Q/ii/ — 3 Jnv — 2K1Zn tJ — TZK.ju + 2K pG TZnuav + / ER-u.\K l 


(2.7) 


( 2 . 8 ) 


where the calligraphic tensor TZ^ pva is the Riemann curvature of 7 ^. For convenience, we 
set the overall scale of the stress tensor to one, l = I 67 tGn = 1 . 

The held equations for the metric g\tN are 


_ 1 1 
0 = Emn — Rmn — ~^9mnR — §9mn — ^QmnLgb 


+ 2 a ( RmabcRn ~ 2RmanbR A — 2 RmaRn + RRmn ) 


(2.9) 


A black hole solution with a hat boundary was found in [58] following previous work [59]. 
In the ingoing Eddington-Finkelstein coordinate, the metric is 


ds z = 2N#dvdr — N^r z f(hr)dv z + r z 5ijdx l dx :1 , i.j = 1, 2, 3. 


To linear order in a, we have 


TV# = 1 — a + 0(a 2 ), 
f(r) = l-\ + 2a (l + 1 ) + 0 (a 2 ). 


( 2 . 10 ) 


( 2 . 11 ) 


Thermodynamics of the EGB black holes was analyzed in [58, 59] . The horizon radius rjj 
and Hawking temperature T are 


1 — a 


r H = 


b 


T = 


1 - 2 a 
7rb 


( 2 . 12 ) 















The conformal boundary is at r = oo. 

To construct fluid dynamics from gravity, we follow [4], First, the static black hole 
geometry (2.10) is boosted along boundary directions x a with a constant boost parameter 
u fl . Then, u ^ and b are promoted into arbitrary local functions of x a , resulting in an 
inhomogeneous geometry 

ds 2 = —2N#u fl (x a )dx fJ 'dr — N#r 2 f (b(x")r) u /1 (x a )u u (x a )dx fl dx u + r 2 V^dx^dx u , (2.13) 

where (x a ) is identified with the fluid velocity and is normalized as ri^ l 'u^( y x a )u 1/ {x OL ) = 
— 1. In general, the metric (2.13) no longer solves the field equations (2.9). Suitable metric 
corrections to (2.13) are needed to make (2.9) satisfied. These corrections are dual to 
parts of II/^jA. Instead of the order-by-order boundary expansion [4], we will collect the 
derivatives in a unified way, as proposed in [14-16] to resum all order linear structures in 
T lw . We linearize u^(x“) and b(x“) 

u^(x a ) = (-l,eui(x a )), b(x“) = b 0 + ebi(x Q ), (2.14) 


where e is an order-counting parameter to be set to unity at the end. Subsequent calcu¬ 
lations are accurate up to linear order in e. The constant bo corresponds to equilibrium 
temperature. For convenience of calculation we set bo = 1. This is equivalent to setting 
7 tT = 1 — 2a, whereas eventually we would like to present our results in units of ttT = 1. 
This is easily achieved by rescaling all the momenta by the corresponding (1 — 2a)-factors. 

The linearized version of (2.13) is 

ds 2 ee d =2N#dvdr — N#r 2 f(r)dv 2 + r 2 5ijdx l dx^ 

—e 1 2 N#Uidx l dr + 4 N# ^1 — ~jdv 2 + 2r 2 [l — N#f(r)] Uidvdx 1 
which is referred to as a seed metric. Formally, we write the full metric as 

ds 2 = gMNdx M dx N = ds 2 ccd + ds^ orr , (2-16) 



where ds 2 orr represents metric corrections. We choose a “background field” gauge [4] 


g rr — 0, g r ^ oc u 


Tr 


,(o) 


-l 


,(i) 


= 0 , 


(2.17) 


where corresponds to the first line in (2.15) and g W denotes metric corrections. Un¬ 
der (2.17), ds 2 orr can be parameterized in the form 

ds 2 0^ = e | dv 2 — 3N#h,dvdr + r 2 hdx 2 + 2r 2 [1 — f(r)] jidx l dv + r 2 aijdx 1 dx^ , (2.18) 

where a tJ is a traceless symmetric tensor of rank two. The functions h, k, ji and aij depend 
on the holographic coordinate r and, through the field equations (2.9), are functionals of 
the fluid velocity u^. 

Boundary conditions for the metric corrections were discussed in details in [15]. The 
first one is that all the metric components in (2.18) are required to be regular over the 
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whole range of r. Second, since the boundary metric is fixed to be r]^, near r = oo we 
demand 

h < O (r°) , k < O (r 4 ) , < O (r 4 ) , a l3 < O (r°) . (2.19) 

Finally, the fluid velocity is defined in Landau frame 

u ,l T IJU = -£ Uv =u»n { ^ = 0 . ( 2 . 20 ) 

Under these boundary conditions, expressions for T^ v greatly simplify. We summarize them 
in Appendix A. 


3 From gravity to fluid dynamics 

In this section, we derive the stress-energy tensor of the boundary fluid by solving the 
field equations (2.9). There are fourteen independent components, which are split into ten 
dynamical equations and four constraints. As in [14-16], our strategy will be to first solve 
the dynamical equations, without imposing the constraints. This turns out to be sufficient 
to uniquely fix the transport coefficients, or in other words we construct an “off-shell” stress- 
energy tensor of the dual fluid. The remaining four constraints are the conservation law of 
the stress-energy tensor. This equivalence is demonstrated in Appendix A. 

3.1 Deriving the fluid dynamics 

The dynamical equation E rr = 0 yields 

(l — 4ct + 4a?’ -4 ) (5d r h + rdf h) = 0. (3.1) 

The asymptotic constraint h < O (r°) and Landau frame convention IL^oo) = 0 lead to 
h = 0. The dynamical equation for k is read off from E rv = 0, 

0 = 3 r 2 d r k — 6 r 4 0u — r 3 d v du + 2dj + rd r dj + r 3 didjaij 

+ ° A [—48 rk + 3 (4r 2 — 3r 6 ) d r k + 8 (3r 8 + r 4 ) du + (5r 7 — 4r 3 ) d v du (3.2) 

—4 (3r 4 + 5) dj — 2 (3r 5 — r ) d r dj — (5 r‘ + 4 r 3 ) d % d 3 a l3 \ , 

which will be solved by direct integration, once solutions for j t and o^j are obtained. 

From E r i = 0, we arrive at the dynamical equation for ji, 

0 = rdfji — 3 d r ji + r 3 d r djOtij + rd 2 Ui — rdidu + 3r 2 d v iii 

- ^ [(5r 5 - 2r) df ji - 3 (5r 4 - 2) 8 r ji + 4 (r 7 + r 3 ) d r d 3 a l3 
+ (5r 5 + r) (d 2 Ui — didu) + 4 (3r 6 + r 2 ) d v Ui] , 
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(3.3) 


which is coupled with ctij only. For the tensor mode aj.j, we find it more convenient to 
consider the combination E t j — ^dijEkk = 0, 


0 = (r 7 - r 3 ) drOtij + (5?’ 6 - r 2 ) d r a^ + 2r 5 d v d r aij + 3 r 4 d v aij + r 3 [[a]]^ 

+ (1 - rd r ) [[j}\ij + (6r 4 + 2r 3 <%) a i3 - ^ [2 (r 11 - 3r 3 ) d 2 a, 

+2 (5r 10 + 9r 2 ) d r a^ + 2 (3r 9 + 4r 5 ) d v d r otij + (9?’ S — 4r 4 ) d v atij ( 3 - 4 ) 

+4 (r 7 - 3r 3 ) [[a]]ij + 5 (r 4 + 6) [\j}]ij - (5r 5 + 6r) \[j}} i3 
+ (18r 8 - 8r 4 ) o^- + 8 (r 7 + r 3 ) ^cr^-] , 

where the notations [[a]]jj, [[j]]ij and <r,y are dehned as 

_+JJ ij = 5 Otjj f didkOtjk + djdkOtik ijdkdlCXklj j 

V 3 / (3.5) 

2 2 K 
[[.)]]■?'j = dijj -F djji y dj• 2<Tij — diUj + djUi -^dijdu. 

Notice that, as in [15], source terms in (3.3,3. 4) are only constructed from m. To solve 
these partial differential equations, we first decompose j t and aij in a basis formed from m, 


ji = a (d v ,d 2 , r) m + b (d v ,d 2 ,r) d t du , 
otij = 2c (<%, d 2 , r ) a %3 +d(d v , d 2 , r) TT tJ , 


(3.6) 


where Oij is dehned in (3.5) and iTij = didjdu—^5ijd 2 du. Then, in Fourier space, dynamical 
equations (3.3,3. 4) translate into a system of second order ordinary differential equations 
for the decomposition coefficients 

0 = rd 2 a — 3d r a — q 2 r 3 d r c — q 2 r — 3iur 2 — [(5r 4 — 2) ( rd 2 a — 3d r a) 

—4 q 2 (r 7 + r 3 ) d r c — q 2 (5r 5 + r) — 4 iG: (3?’ 6 + ?' 2 )] , 

0 = rd 2 b — 3d r b — -q 2 r 3 d r d + -r 3 d r c — r — \ (5?’ 4 — 2) (rd 2 b — 3 d r b) 

3 3 r 4 

— - (r 7 + r 3 ) d r (2 q 2 d — c) — (5r 5 + r) , 

3 

0 = (r 7 — r 3 ) d 2 c + (5r 6 — r 2 ) <9 r c — 2iur 5 d r c — rd r a + a — 3iu)r 4 c 

— iCur 3 + 3r 4 — [2 (r 11 — 3r 3 ) <9 2 c + 2 (5r 10 + 9?’ 2 ) d r c 

—2icD (3r 9 + 4r 5 ) <9 r c — iui (9r 8 — 4r 4 ) c + 5 (r 4 + 6) a 

— (5r 5 + 6r) <9 r a + (9r 8 — 4r 4 ) — 4 iCo (r 7 + r 3 )] , 

0 = (r 7 — r 3 ) <9 2 d + (5r 6 — r 2 ) <9 r d — 2 iur 5 d r d — -^r 3 (2c — q 2 d) + 2b 

o 

— 2rd r b — 3iur A d — ° 4 [2 (r 11 — 3r 3 ) d 2 d + 2 (5r 10 + 9r 2 ) d r d 

—2iCd (3 r 9 + 4r 5 ) d r d — iCo (9r 8 — 4r 4 ) d — ^ (r 7 — 3r 3 ) (2c — g 2 d) 

+10 (r 4 + 6) 6 — 2 (5r 5 + 6r) d r b\ . 
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Equation (3.2) becomes 


( 4a \ , 16a, ( 2 1 _ 2 , _ 2 . 

0 = [l - 3a + <%£: - -^-A; - |2r - -zwr - — ( a - q b) 


-fr ( °r a - fdrb) - -<i 2 r (fd -2 c) + ^ 


9 


a 


-*w(5r 5 — 4?’) 


+ (3r 4 - l) ( d r a - q 2 d r b ) + ^-q 2 r (5r 5 + 4r) (g 2 d - 2c) 

or y 

+ ^2 ( 3r4 + 5 ) ( a - “ jj ( 3r6 + r2 ) } 


(3.8) 


The barred momenta are defined as Q = (1 — 2a)cc and q = (1 — 2a)q, which emerge as a 
result of the above mentioned rescaling of units. 

We first study the large r behavior of the metric corrections, which propagates into 
the expression for the fluid’s stress tensor. The velocity dependence of T^ v enters via the 
decomposition (3.6). Examining equations (3.7) near the conformal boundary r = oo, it is 
straightforward to show that 


a ■ —iu (1 + a)r 3 + O , 
cJ ^oo > 1 -a + C*(u q 2 ) + p/l\ 

y'* \ jit) J 


d r^ D ${u,q 2 ) +o 



(3.9) 


where C * 4 and D£ are unknown coefficients, which cannot be determined from the asymptotic 
analysis alone. To compute them, we have to integrate (3.7) over the entire bulk. Regularity 
of the metric components in (2.18) imposes two boundary conditions at r = rn, which are 
sufficient to fix and D 4 uniquely. The large r behavior of k is 

(1 — 2a) iter 2 j> du + O 

Boundary conditions (2.19,2.20) were imposed in deriving (3.9,3.10). 

Plugging (3.9,3.10) into (A.l, A.2, A. 3), we obtain the boundary stress-energy tensor 


(() • (3.10) 


, ,2 , N o 2 

-(!-a)r 3 + - 


' T 00 = 3 (1 - 5a) (1 - 4bi), 

T 0 i = Ti o = -4(1- 5a) m , 

Tij = 5ij (1 - 5a) (1 - 4bi) 

+ 4(1- 3a) [2 C 4 (u), q 2 ) (1 + 6a) + T> 4 (a), q 2 ) (1 + 2a) vr^] . 


(3.11) 


Covariantization of (3.11) gives standard expressions (1.2,1.6) of T^, with e and P given 
by (1.12). The viscosity functions q and £ re-expressed in units of irT = 1 are 


77 (w, q 2 ) = —4 (1 — 3a) (7 4 [(1 — 2a) w, (1 — 4a) q 2 ] (1 + 6 a), 
c (u, q 2 ) =-4(1- 3a) D$ [(1 - 2a) w, (1 - 4a) q 2 ] (1 + 2a). 
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3.2 Gauss-Bonnet corrections to the viscosity functions 

To determine the viscosity functions, we are now ready to fully solve the dynamical equa¬ 
tions (3.7). In the next subsection, we start with the hydrodynamic limit and solve (3.7) 
perturbatively in momenta. In this way, we reproduce some known results in the literature 
and also obtain a set of new third order transport coefficients. In the subsection to follow, 
we address our main goal of resumming all-order derivative terms. This will be achieved 
by numerically solving (3.7) for generic values of a) and q 2 . 


3.2.1 Analytical results: hydrodynamic expansion 

We introduce power counting parameter A by w -> Aw and q^ —> Xqi, and expand the 
decomposition coefficients (3.6) in powers of A, 


OO 

a (w, qi,r) = 22 A n a n (w, r), 

n =0 
oo 

c (w, qt, r) = 22 ^ c n (w, qi , r), 

n =0 


oo 

b (w, qi, r) = 22 (<+ qi, r ), 

77 .— 0 
OO 

d (w, qi, r) = 22 ^ n dn (<+ qi, r ). 

77.—0 


(3.13) 


At each order in A, there is a system of ordinary differential equations for a n etc, whose 
solutions are double integrals. In Appendix A, we summarize these results. Then, C^ and 
-Dl) are expanded as 


C fe 4 (w, q 2 ) = - - (1 - 8 a) - Ucb [(2 - In 2) - (23 - 6 In 2) a] + q 
4 8 


i 


" 2 - 0.497(227)a 

OZi 


+ w 2 {y^ ( 67 t — 7 r 2 + 24 — 36 In 2 + 12 In 2 2) - 1.56(140) a| + ••• , (3.14) 


D^(uj, q 2 ) =— (7r — 5 + 2 In2) + — (41 — 7tt — 14In2) a + • • • , 


1 


1 


48 


24 


where in C 4 we have only numerical results for the linear in a second order terms. The 
viscosities (1-14) are obtained by substituting (3.14) in (3.12). 

Taking plane wave ansatz for m and bi, the conservation law d IJ 'T fIl/ = 0 results in 
dispersion equations 


shear wave : 
sound wave: 


(1 + 3a) w + Uq 2 T] (w, q 2 ) = 0, 

(1 + 3a) (g 2 — 3w 2 ) — iujq 2 r] (w, q 2 ) + -iujq^C (w, q 2 ) = 0. 


(3.15) 


In the hydrodynamic limit, the dispersion equations (3.15) could be solved perturbatively. 
For the lowest modes they read 


shear wave : 
sound wave : 


w =- | (1 -8a)q 2 - ^ [1 -log2+ (8 #i - 40 + 16log2) a] q 4 H-, 

u = ± -2= - - (1 - 8a) q 2 ± -?-= [3 - 2 log 2 + (16 log 2 - 38) a] q 3 

VS 6 V 24\/3 ^ f (3.16) 

- k 2 - 24 + 24 lQ g 2 - 12 lQ g 2 2 + (#1 + 444 #2 - 294 
864 

- 90vr - 3vr 2 + 60 log 2 + 36 log 2 2) a] q A + ■ ■ ■ , 
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where #q = 6.53(280) and ^2 = 2.11(320) are known numerically only. These hydrody¬ 
namic modes should agree with the lowest quasi-normal modes of the EGB gravity. 

3.2.2 Numerical results: all-order resummed hydrodynamics 

For generic values of oj and q 2 , we resort to a shooting technique and solve (3.7) numerically. 
Our numerical procedure is essentially the same as that of [15]. We start with a guess 
solution at the horizon r = th and integrate (3.7) until the conformal boundary r = 00 . 
Then, we fine-tune the initial guess until thus generated solution satisfies the boundary 
conditions at r = 00 . 



Figure 2. The viscosity 77 as function of w and q 2 . 


Numerical results for the viscosities are shown as 3D plots in Figures 2 and 3, and then 
sliced at q = 0 or u = 0 in Figures 4, 5 and 6. A marking behavior of all the functions is 
that they vanish at very large momenta, a behavior necessary for restoration of causality. 
Damped oscillations are clearly visible reflecting a complex pole structure of the viscosities 
as functions of complex uj. These are the quasi-normal modes of the so-called scalar (or 
tensor) channel [60-62], 

Without the Gauss-Bonnet corrections, the viscosities and display only a weak 
dependence on spatial momentum q. meaning the dissipation is quasi-local in space. In 
contrast, j/ 1 ) and introduce a much more noticeable space dependence. In order to 
see relative correction to viscosity function due to Gauss-Bonnet term, in Figure 7 we 
combined r/ 0 - 1 and r/ 1 ) for upper and lower bounds of a. We observed that the Gauss-Bonnet 
correction does introduce a profound spatial dependence for the viscosity. In addition, 


- 12 



Re[f (1) ] 



«m[£ (1) ] 


Figure 3. The viscosity ( as function of ui and q 2 . 


Re[^ (0) ] 
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Figure 4. The viscosity 77 as function of u; with q = 0. 


the Gauss-Bonnet correction changes shape of the viscosity noticeably in the intermediate 
regime of momenta where the amplitude of oscillation is enhanced. 

Another interesting observation concerns imaginary parts of rj. While Im[ry(°^] is always 
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Figure 5. The viscosity £ as function of u> with q = 0. 



Figure 6. The viscosities rj and £ as functions of q 2 with u> = 0. 


positive, Imf?/ 1 )] changes sign. This implies that for certain values of a, both positive and 
negative, Im[r;] may become negative. To clearly see this behavior, in Figure 8 we plot 
Im[r/] as function of uj (q = 0 ) and a (within the causality interval ( 1 . 11 )). Im[^] becomes 
negative when a goes below the critical value —0.05. With q increased, this critical value 
gets larger. If the viscosity function had an interpretation of a correlation function, then its 
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Figure 7. The viscosity rj as function of w and q 2 for a = 9/200 (top) and a = —7/72 (down). 



Figure 8. The viscosity rj as function of u (q = 0) and a within the causality interval (1.11). The 
hole around oj = 4 indicates the region where Im[?y] < 0. 

imaginary part would be a spectral function and would have to be positive. Yet, beyond the 
first order in the gradient expansion the correlation functions get additional contributions 
from so-called gravitational susceptibilities of the fluid [3, 16, 17]. So, while the possibility 
that Irn[r;] becomes negative for some values of a does not immediately imply a problem, 
we take it as a signal for possible issues with causality in the theory. 
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To better explore the effect of the Gauss-Bonnet corrections, we now represent our 
results as memory functions in real time. Let perform an inverse Fourier transform of 
r](uj,q 2 ) with respect to co only, 

/ °o J. , 

(3.17) 

-oo V 



Figure 9. Memory function fi{t,q 2 ). Left: 3D plot as a function of time t and momentum squared 
q 2 . Right: 2D plots as function of time t: different curves correspond to different q 2 ( q 2 = 0,1, 2, 3 
from the rightmost for and from the bottommost for 

In Figure 9, we plot the time dependence of the memory function fj(t, q 2 ). As has been 
pointed out in Introduction, fj(°> vanishes for negative times, consistency with causality 
requirement. The Gauss-Bonnet correction also has support in positive times only. 
Similar effect is also found for the second memory function (. 

4 Summary and discussion 

In this work, we discussed effects of Gauss-Bonnet corrections on holographically dual fluid 
dynamics. For bulk EGB gravity, we found a boosted black brane solution with locally 
perturbed horizon. Our construction is accurate to linear order in amplitudes of fluid 
velocity temperature T(x), and the Gauss-Bonnet coupling a. This black brane 

solution is dual to all-order linearly resummed fluid dynamics, and a-corrected viscosity 
functions were read off from it. 
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In the hydrodynamic limit, we reproduced known results for a-corrected shear viscosity 
and relaxation time. As a new result, we expanded the knowledge on transport coefficients 
by computing a corrections to the third order coefficients. Beyond the hydrodynamic limit, 
we computed Gauss-Bonnet-corrected viscosity functions. We observe two qualitatively 
new effects induced by the corrections. First, the viscosities become less local in real space. 
Second, due to a corrections, Imf^j can become negative. 

Finally, we Fourier transformed the viscosity functions into real time, where they play 
a role of memory functions. For positive times, we observed a pattern of damped oscillation 
reflecting a structure of complex poles. Interestingly, the poles of are apparently shifted 
compared to the ones of rj^\ Given that the EGB memory function does not display any 
causality violation, we do not expect any dramatic a-induced effects on EGB quasi-normal 
modes. However, we think a study of EGB quasi-normal modes might provide additional 
insight on the problem. This is, however, beyond the scope of the current paper. 


A Computational details 

In this Appendix, we provide some computational details which were omitted in deriving 
the boundary fluid dynamics. 

In terms of the metric corrections (2.18), the tensor T^ v is 
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Tij =5ij (1 — 5a) (1 — 4ebi) + eSij r 3 du + 9 r A h + 2 r 5 d r h — r 3 d v h + -r 2 d 2 h 
—r 2 d 2 h + k — rd r k + -d v k -3—<9 2 k — -dj + \d v d j — -r 2 d k dia k i 


+ e 


2 r 2 

T^didjk - r 3 (diUj + djUi) - ^ r 2 didjh - r 5 d r a^ ^r 2 d 2 otij 


v d v (Xij “I - ^ ^ 2r 2 ^0 djji) T 2 r ~ (did k Q!j k T djdk&ik) 


+ o r dv<*ij 


— eaSi-j ( 8 r 3 du + 3 k — 3rd r k H—<9„&-— 27r A h 

\ r 2 r 2 

+3r 5 d r h + Ar 3 d v h + \r 2 d 2 h — br 2 d 2 h + ~^d v d j — —dj — ^ r 2 d k dia k i 


(A.3) 


— ea 


2 v r 2 “~ J r 2 

T^didjk - 4r 3 (djUj + dfu,) - ^ r 2 didjh + - ^9,^ {dijj + djji ) 


5 5 5 

3r d r otij 4r d v otij — v d oiij + — r d v Oiij -|- — r (did k otj k T djdkomj) 


where we have dropped terms that explicitly vanish at r = oo. 

For consistency, constraints in (2.9) have to be satisfied by the gravity solution pre¬ 
sented in section 3. We find it more convenient to consider suitable combinations of 
Emn = 0. The first one is E vv + r 2 f(r)E vr = 0, 


0 = 4r 3 9tt + r 2 d 2 h\ — \2r 3 d v h\ — ( r 2 + r 6 ) d r du — 4r 3 9j — r 2 d 2 k + 3 r 3 d v k 

+ 2r 2 d v dj + (r 4 — l) d r dj — [2r 5 (3r 4 + 7 ) du + 24r 4 9 2 bi — 60r 5 9<,bi ^ 

— (3r 8 + 3r 4 — 2) d v du — 2 (Hr 5 — r) dj — 4 (r 4 — l) d 2 k — (r 8 — r 4 ) didjCtij 
+3 (3r 5 — 4r) d v k + 4 (3r 4 — l) d v dj — 3 (r 7 — r 3 ) d r k + 4 (r 6 — 2r 2 ) 9 r 9j] . 

The combination E v i + r 2 f(r)E r i = 0 yields 

0 = r A d 2 Ui — r A didu + 4r9jbi — 4 rd v Ui — r A d 2 Ui — d 2 ji + d % dj — r9*fc + 4 rd v ji 
+ r A d v d k a ik + (r 6 - r 2 ) d r d k a ik + r 2 d r dik - r 2 d r d v ji - ^ [(3r 9 - 23r 5 ) 9,/(q 

+20r 5 9jbi + (5r 8 + 3r 4 ) ( d 2 Ui — d{du) —4 (r 8 — r 4 ) 9 2 ttj — (3r 5 — 20r) d t k ^ ^ 
-6 (r 4 + l) (9 2 j* - d,dj) + 20 r 5 d v j t - 3 (r 7 - r 3 ) d r j t + (3r 6 - 4r 2 ) 9 r 9,A; 

- (5r 6 - 2r 2 ) d r d v ji + (r 8 - r 4 ) 9 2 j; + 2 (2r 10 - r 6 - 3r 2 ) d r d k a ik 
+4 (r 8 + r 4 ) 9„9 fc a« : ] . 


With the near r = oo behaviors (3.9,3.10) at hand, the large r limit of (A.4,A. 5) can be 
shown to produce the conservation law d^Tu v = 0. 

In the hydrodynamic limit, we perturbatively solved holographic RG flow equations (3.7). 
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Recall the formal expansion (3.13) 


a (uj, qt,r ) = ^ ( w > Qu r ) > & (w, Qi, r ) = ( w ’ Qi ’ r ) ’ 


n=0 

oo 


n=0 

oo 


(A. 6 ) 


c (w, qi, r) = ^2 ^ n Cn (w, qi, r) , d (uj, qi, r) = ^ A"d n (w, qi, r ). 

n=0 n=0 

Then, perturbative solutions for the metric corrections can be expressed as double integrals. 
We here summarize the main results, 


dx 


co = 


J r (x 5 — x) — 2a (x 5 — 3x -3 ) 

r—»oo 1 ct 1 8 ck „ / 1 \ 

--4^ + °U ’ 


a 0 = 0, (A.7) 

[-3 y 2 + a (9 y 2 - 4 y~ 2 )] dy 


' r H 


a i = — ico 


2a 


ci = 


i; 


dx 


(1 + a) r 3 H- 


dy {3 iujy 2 c 0 (y) + 2iuy 3 d y c 0 (y) - iuy 


(A. 8 ) 


(A.9) 


x 5 — x — 2a (x 5 — 3x -3 ) JrH 
-aid: [(9 y 2 - 4 y~ 2 ) c 0 (y) + 2 (3 y 3 + 4 y _1 ) d y c 0 (y ) - 4 y + 12 y~ 3 ] } (A.10) 


-^[ 2 -l" 2 -«(23-61n2)]+o(4 


r /»oo 

3 


bo = — x 6 dx 


dy 


y- y 3 d y c 0 (y)/3 + a[4(y 3 + y 3 ) d y c 0 (y) - (by + y 3 )] 


Jrn J x 

3 2 r-s>oo 1 2 , /n / - i 

-8 + 3 a — + °U. 1 ' 


y A -a (5 y 4 - 2) 


(A.ll) 


/■°° , r 

a 2 = / dxx 

Jr J x 


3 / d ^ q 2 y 3 d y co(y) + q 2 y- a [4q 2 (y 3 + y l ) d y c 0 (y) + q 2 (by + y 3 )] 


y 4 - a (5y 4 - 2) 


r^J_ }2(1 _ 7a) + 0 (J_ 


do — — 


[ 


*- x - 2a ; 5 _ ss-3) /„ + ^ c » fe) 


—a 


y 

10 12 \ _ , , , /10 60 \ , . 8 , 0 _ 3 . , . 

7 + y5 J d v b o(y) -\y2 + tf}J &o(y) + 3 O' - ) c °(y) 


~ 4^4 [6 — tt — 2 in 2 — a (82 — I4ir — 28 in 2)] + O f j , 


C2 = - 


i: 


dx 


(A.12) 


(A.13) 


dy {2iuy 3 d y ci(y) + 3iuy 2 ci(y) + y 1 d y a 2 (y) 


x 5 — x — 2a (x 5 — 3x 3 ) Jr 
~y~ 2 a 2 (y) ~ a [2iw (3y 3 + 4y _1 ) <9 y ci(y) + id (9y 2 - 4y“ 2 ) ci(y) ''' ' ‘' 
-b (y~ 2 + 6 y~ 6 ) a 2 (y) + (5y _1 + 6 y~ 5 ) 9 y a 2 (y)] } , 
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where rjj = 1 — a as defined in (2.12). From large r behavior of these functions, we arrive 
at the power expansion (1.14) of the viscosity functions. 
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